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ABSTRACT 

Aims. We attempt to determine the nature of the high energy emission of the radio galaxy 3C 111, by distinguishing between the 
effects of the thermal and non-thermal processes. 

Methods. We study the X-ray spectrum of 3C 111 between 0.4 keV and 200 keV, and its spectral energy distribution, using data from 
the Suzaku satellite combined with INTEGRAL, Swift/BAT data, and Fermi/LAT data. We then model the overall spectral energy 
distribution by including radio and infrared data. 

Results. The combined Suzaku, Swift and INTEGRAL data are represented by an absorbed exponentially cut-off power-law with 
reflection from neutral material with a photon index T = 1.68 ± 0.03, a high-energy cut-off E cut = 227^^ 3 keV, a reflection component 
with R = 0.7 ± 0.3 and a Gaussian component to account for the iron emission-line at 6.4 keV with an equivalent width of E W = 
85 ±11 eV. The X-ray spectrum appears dominated by thermal, Seyfert-like processes, but there are also indications of non-thermal 
processes. The radio to y-ray spectral energy distribution can be fit with a single-zone synchrotron- self Compton model, with no need 
for an additional thermal component. 

Conclusions. We suggest a hybrid scenario to explain the broad-band emission, including a thermal component (iron line, reflection) 
that dominates in the X-ray regime and a non-thermal one to explain the spectral energy distribution. 
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1. Introduction 

Radio galaxies are a subclass of active galactic nuclei (AGNs) 
and in the unification model of AGN radio galaxies are the 
radio-loud counterpart s of Seyfert galaxies (lAntonuccil 1 1 993b 
lUrrv & Padovanilll995l) . For both of these classes, it is thought 
that the inclination is large, such that the observer views the core, 
which is dominated by thermal processes through absorbing ma- 
terial. This differents from blazars, where the inclination angle is 
very small. The observed emission originates in the relativistic 
jet, which is dominated by non-thermal processes. Owing to the 
Doppler boosting of the emission in the jet, the energy of the ob- 
served emission from blazars can reach the y-ray regime, even 
up to TeV energies. Several non-blazar AGNs have also been 
found to emit significantly in the y-ray regime (IHar tman et al. 
2008; lAbdo et all 1201 Oct lAckermann eTatl |20111). Although 
there have been several theories (i.e. misaligned jet, shocks in 
the radio lobes) about the origin of this high-energy radiation, 
conclusive support of any of t hese theories is yet to be found. 
3C 1 1 1 is a ne arby (z=0.049.ISargendl 19771) flat-spectrum FR- 
II radio galaxy dFanaroff & Rilevll 19741: ILinfield & Perlevll 19841) 
that d isplays sttong, broad emission-lines in the optical (Sargent 



1977) and an iron-emission line in the X-ray regime (ILewis et al.1 
2005), similar to Seyfert galaxies. In the radio regime, the 



galaxy's emission is dominated by the relat ivistic jet, which has 
an angle of 18° to our line of si ght dJorstad et al.l2 005). The pro- 
jected size of the jet is 78 kpc dBridle & Perlevi ll984). There is 
no visible counter-jet, although a bright lobe is detectable in the 
opposite direction of the observed jet, which is likely fed by the 
undetected counter-jet (ILinfield & Perlevll 19841) . 
Earlier studies of 3C 111 in the X-ray band h a ve ide ntified 
the high energy cut-off of the spectrum. Dadina] (120071) found 



a lower limit of £ C ut ^ 82keV using BeppoSAX data and 



IB alio et alj d201 ll) derived a similar lower limit to the cut-off of 
E alt > 75 keV, using data fro m XMM-Newton and SuzakufXIS 
and PIN. Moli na et all (120091) . using XMM-Newton, Swift, and 
INTEGRAL data, constrained the cut-off energy to be E cut = 
126^Q 3 keV. It is uncertain whether a reflection component is 
present in the X-ray spectru m. An upper limit of R < 2.25 
was found by iDadinal d2007l) using Bepp oSAX data and con- 
straints of R = 0.9 + 0.6 were derived b y lMolina et all (120091) 
and R = 35 +0.06 by IB alio et al.1 d201ll) using Suzaku/XIS and 
PIN data. iBallo et al.1 d201 lb used XMM-Newton data to find the 
reflection in the en ergy range 0.4-10 keV to be R = 0.19^^. 
River s et al ] d2011h did not detect a reflection component, us- 
ing RXTE data between 3 keV and > 100 keV. The equiva- 
lent width (EW) of the iron line at 6. 1 keV has been mea- 
sured several t i mes a nd is likely variable. lBallo et all (1201 ll) and 
Tomb esi et al.l (1201 Ol) found similar values of EW — 75 + 13 eV 
{Suzaku/XIS and PIN) and EW = 86 + 16 eV (Suzaku/XIS), re- 
spectively. Using Suz aku/XIS data of three d ifferent observations 
taken two years later. iTombesi et all (1201 ll) found a s maller EW 
of betw een EW > 33 eV and EW > 40 eV. Similarlv. lBallo et al.1 
d201ll) used XMM-Newton data to fit the 0.4-10 keV spectrum 
and found an EW = 38^' eV. 

3C 111 has been detected in y-ray s by CGRO /EGRET 
(IHartman et al.lll999t ISguera et al.ll2005t IHartman et all 12 008). 
and wa s included in the first Fermi/LAT catalogue dAbdo et al.l 
1201 Oal) . In the second Fermi/LAT catalogue, the source was 
omitted since it had no longer been significantly detected 
(Nol an et al.ll2012l) . 3C 1 1 1 is likely variable in the y-ray regime 
(Ackermann et aDl201 ll) . Using 24 months of Fermi/LAT data, 
iGrandi et al.l(l2012l) found 3C 1 1 1 to be detectable during a short 
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time-period ( At ~ 30 - 60 days), limiting the radius of the emis- 
sion region to be R < 0. 1 pc (assuming a Doppler factor of 6 = 3) 
based on to causality arguments. The detectability of the source 
in the GeV energy range coincided with an increase in the flux in 
the millimeter, optical and X-rays regimes, indicating the emis- 
sion is likely to emerge from the same region. Since the outburst 
in millimeter, optical, and X-rays is associated with the ejection 
of a bright radio knot, this indicates that the GeV emission orig- 
inates from the radio core within 0.3 pc of the central supermas- 
sive black hole. 

To understand the physical processes causing the emission in 
3C 111, we first study the X-ray to y-ray spectrum to find 
whether the emission is the product of thermal (Seyfert-like) or 
non-thermal processes (blazar-like) or a combination of both. 
We then model the spectral energy distribution (SED) to un- 
derstand the processes dominating the broad-band radiation of 
3C 111 from radio to high energies. 



2. Observations and data reduction 

To cover a large energy band from X-rays to y-rays and study the 
high-energy emission from 3C 111, we used data from several 
different instruments. 



2.1. Suzaku-XIS/PIN 

3C 111 was observed by Suzaku from 22 August 2008 to 25 
August 2008 with a total elapsed time of 236.9 ks in HXD 
nominal pointing mode. W e analysed data fro m the X-ray 
Imaging Spectrome ter (XIS, Koyama et al. 2007) and Hard X- 
ray Detector (HXD, Takahas hiet al.ll2007l) . The XIS instrument 
consists of three separate CCD detectors and has an energy range 
of ~ 0.2 - 12.0 keV. The HXD is a collimated instrument, con- 
sisting of two independent detector systems, silicon PIN diodes 
that function in the range ~ 10-60 keV and GSO scintillation 
counters covering the energy range ~ 30 - 600 keV. 
For the observations made with the XIS, we used the clean 
events provided by the instrument team where the standard event 
cuts had been applied. 

The flux of 3C 111 is too low above 70 keV to extract a sig- 
nificant spectrum from the HXD/GSO detector, hence we used 
only HXD/PIN data. The Suzaku team provides a response file 
for the PIN, which depends on both the epoch in which the data 
were taken and the pointing mode. In addition a non X-ray back- 
ground (NXB) per observation is available for the PIN analy- 
sis. After filtering the source and background files in time using 
the good time intervals, we extracted the source and background 
spectra. The source spectrum needed to be corrected for dead- 
time events, which we achieved using hxddtcor. The NXB was 
simulated with ten times more counts to suppress statistical er- 
rors. We therefore needed to increase the exposure time of the 
background spectrum by a factor of ten. Using xspec, we sim- 
ulated the cosmic X-ray background (CXB) using the PIN re- 
sponse file for the flat emission distribution in the proper epoch 
and point mode. We combined the NXB and CXB files and used 
this as a background in our analysis. 



2.2. INTEGRAL IBIS/ISGRI 

The Internati onal Gamma-Ray Astrophysics Laboratory 
(INTEGRAL; Win kler et all (120031) 1 is a y-ray observatory with 
several instruments on board. We used the INTEGRAL Soft 
Gamma-Ray Imager (ISGRI), which is part of the Imager on 



Board INTEGRAL Spacecraft (IBIS), a coded- mask detector 
ISGRI is sensitive between 15 keV to 1 MeV dLebrun et all 
l2003h . 

We used all data collected by INTEGRAL since the launch of 
the satellite up to August 2009, with a total exposure time of 508 
ks. We first created individual spectra for all science windows 
(a science window is all data produced during one pointing) 
using the ibis_science_analysis routine. We then summed all the 
individual spectra to achieve a higher signal-to-noise ratio. 

2.3. Swift-BAT 

The Burst Alert Telescope (BAT) aboard the Swift satellite is a 
coded-ap erture camera with an energy range of 14-195 keV for 
imaging (Barthelmv et al. 2005). Since the BAT telescope mon- 
itors the sky continuously and has a large field of view (1.4 sr, 
partially coded), it has regularly observed 3C 111. 
We used data from the 5 8 -month hard X-ra y survey, from 2004 
December to the end o f 2009 May (see iTueller et al.1 1201Q- 
Baumgart ner et alj|2010l) . The detected flux for 3C 111 in this 
survey is 1.2 x 10 _10 ergcm _2 s _1 between 14keV and 195 keV. 



2.4. Fermi -LAT 

The Large Area Telescope (LAT) aboard the Fermi satellite 
operates in an energ y range between 20 MeV and 300 GeV 
(lAtwood et al.1 12009 ). The LAT is a pair-conversion telescope 
with a very wide field of view that scans the sky continu- 
ously. During the nominal all-sky survey observing mode of 
FermifLAT, a total effective exposure time of 83.7 Ms was ac- 
cumulated in the direction of 3C 111 between 4 August 2008 
and 20 April 2011. 

We selected diffuse-event-class photons (P6_V3 instrument re- 
sponse functions) between 100 MeV and 200 GeV in a cir- 
cular region of interest with radius of 15° around the source. 
Events with a zen ith angle of more than 105° were excluded 
dAbdo et alj|2009l) . We used the standard cuts proposed by the 
Fermi team based on the data quality of the events and the in- 
strument configuration. 

The maximum-likelihood analysis tool gtlike models all the 
emission in a given region, which can contain several sources. 
The goodness of fit is expressed as the log-likelihood value 
which is the probability of obtaining the data given an input 
model. After creating a model including all detectable sources in 
the field, we used gtlike to fit the model until the log-likelihood 
value was maximised. We found the test statistic value of 3C 111 
to be TS = 12.8, which corresponds to a significance of ~ 3cr. 
The data provided only two significant points so it is impossi- 
ble to fit the spectrum in e.g. xspec. However, when creating the 



spectrum, gtlike gives a spectral slope of F y = 2.41 ±0.17 across 
the range 0.1-200 GeV and a flux off = 1.2 x 10~ 8 phcirT 2 s _1 . 



3. Results 

3.1. X-ray spectrum 

We simultaneously fitte d the X-ray sp ectra between 0.4 keV 
and 200 keV, in xspec (lArnaudl [l9961 We started our analy- 
sis by fitting the spectrum with an absorbed power-law and a 
Gaussian component at 6.1 keV to account for the redshifted 
iron emission-line. The fit resulted in a reduced chi-squared of 
xl - 1.14 (1678 d.o.f.). We tried to improve our fit by using 
a cut-off power-law, which yielded a better fit with xl — 1-12 
(1677 d.o.f.). We explored the possibility of a broken power-law, 
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but with axl — 1.14 (1676 d.o.f.), this gave no improvement in 
the fit. Adding a reflection component did further improve the fit, 
X 2 V = 1.10 (1676 d.o.f., F-test probability 4 x 1(T 8 ). The result- 
ing fit and the residuals are shown in Figure [TJ We derived the 
following best-fit parameters and 90% errors: the value for the 
equivalent hydrogen column-density A^h = (9.0 + 0.2) x 10 21 
cm -2 , a power-law index of F = 1.68 + 0.03, a high-energy 
cut-off E cat = 227^g!j 3 keV and a reflection scaling factor of 
R - 0.7 + 0.3. Figure [2] shows the confidence contours of the 
reflection component versus the cut-off energy with Ay 2 of 68%, 
90%, and 99.7%. We can exclude a cut-off below 130 keV at the 
99.7% confidence level. 

We fitted the data with a more physical model, compPS 
dPoutanen & Svenssonlll996h . which describes the process of 
thermal Comptonization. In this model, seed photons from the 
cool, thick accretion disc are injected into the electron plasma. 
The electron cloud can have several geometries, and we chose a 
plane-parallel slab: we also applied other geometries but found 
that they neither influenced the parameter values nor improved 
the fit. The electrons in the cloud have a Maxwellian distribution 
with an electron temperature T e and an optical depth t related 
to the Compton parameter y = 4-TkT e /(m e c 2 ). The seed photons 
are up-scattered from their initial energy E-, to Ef — e y E t . The 
Compton scattering of the seed photons with the electrons in the 
plasma results in a spectrum that is afterwards reflected from a 
cool medium and then smeared out by the rotation of the disc. 
The best-fit model gave a. xl - 1.10 (1676 d.o.f.). The values 
that we found for this model are an equivalent hydrogen column- 
density A^h = (9.0 + 0.2) x 10 21 crrT 2 , a temperature of the elec- 
trons kT e = 91^gkeV, a reflection component R - 1.8*g5, and 
a Compton parameter of y = 0.6 + 0.1, which corresponds to an 
optical depth t = 0.8. 

We also fitted all data sets individually with an absorbed power- 
law and a Gaussian component. The values of the flux levels and 
power-law indices for these individual best fits to the data can be 
found in Table [TJ 

The measurements of the iron line at 6.4 keV are consistent in 
all three Suzaku/XIS spectra. After adding the three spectra to- 
gether, we derived a line energy £);„ = 6.1 1 + 0.02 keV and an 
£W = 85 + lleV. 

3.2. Fitting of the time-averaged SED 

After we had created a SED based on our data, we modelled it 
wit h a public synchrot r on sel f-Compton (SSC) code developed 
by iKrawczvnski et all £20040 The SSC mechanism assumes 
an isotropic population of high-energy electrons that emit syn- 
chrotron radiation followed by inverse Com pton scattering of th e 
synchrotron photons to higher energies dMaraschi et al.N l992). 
The inverse Compton component may also include a thermal 
component of external seed photons (e.g. from the broad-line re- 
gion or accretion disc) that are also Com pton upscattered to high 
energ ies (external Compton component, Dermer & Schlickeiser 
1993). The electron population is located in a spherical volume 
of radius R s with a randomly orientated magnetic field B that 
moves relativistically towards the observer with a bulk Lorentz 
factor of T and an angle between the jet and the line of sight 
of 6. Thus, the radiation is Doppler-shifted with a Doppler fac- 
tor 6 = [r(l - BcosO)]^ 1 . The electron energy spectrum in the 
jet-frame follows a broken power-law, with indices p\ and p2, 
and is characterised by a minimum and maximum energy (E m i n , 
Emax) and a break energy (£\,r)- The power in the SED is mostly 



1 http ://j elley. wustl . edu/multi wave/spectrum/?code 




Fig. 1. The count spectrum of the combined Suzaku/XIS (0.4- 
10 keV), Suzaku/PIN (12-60 keV), INTEGRAL ISGRI (20-200 
keV), and Swift/BAT (15-150) data with the fitted model: an ab- 
sorbed cut-off power-law with reflection from neutral material 
and a Gaussian component to account for the iron line at 6.4 
keV. The bottom panel shows the residuals in terms of the stan- 
dard deviation with error bars of size one sigma. 

500 f ' ' ' 7-H ' ' ' ' 1 ' ' ' fyi : 
400 - / / % * -_ 




100 r 



0.0 0.5 1.0 1.5 

reLrefl 

Fig. 2. Error contours for the best-fit pexrav model, reflection 
component R vs the high-energy cut-off £ cut . The contour levels 
correspond to Ay 2 of the 68%, 90%, and 99.7% statistical confi- 
dence regions. A cut-off below 130 keV is excluded at a 99.7% 
level. 

influenced by the Doppler factor and both the radius and mag- 
netic field of the emitting region. The shapes of the synchrotron 
and inverse Compton peaks and the peak frequencies depend 
strongly on the energies that define the electron population. 
Since we only had data in the X-ray and y-ray regimes, it was 
impossible to model the entire SED. We therefore added archival 
radio data, neglecting the observations that only included the 
core, and infrared data from the NED archival. Assuming that 
the infrared emission is produced in the same region as the X- 
rays (the accretion disc or inner jet), we were able to use the col- 
umn density found in the X-ray spectra, A^h = 9.0 x 10 21 cm -2 
to calculat e the extinction in the V -band Ay = Nn/(1.79 X 
10 21 cm" 2 ) dPredehl &Schmitll 1995b . From the extinction in the 
V-band, we calculated the extinction in th e near-infrared bands J , 
H, and K using the correction factors from ISchlegeletal.l (fl998). 
We removed the Suzaku/XIS data below 1 keV to avoid possible 

2 http://ned.ipac.caltech.edu/ 
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Table 1. The parameters for the individual power-law fits of the data, and their 90% confidence levels. 



Instrument 


Epoch 


Exposure time 


PL index 


/[erg cm 2 s '] 


Energy range 


XISO 


22-25 August 2008 


95.4 ks 


1.60 ± 0.02 


(2.20 ±0.02) x 10-" 


0.4-10 keV 


XIS1 


22-25 August 2008 


95.4 ks 


1.59 ±0.02 


(2.24^-g) x 10- 11 


0.4-10 keV 


XIS3 


22-25 August 2008 


95.4 ks 


1.63 ±0.02 


(2.2 ± 0.02) x 10-" 


0.4-10 keV 


PIN 


22-25 August 2008 


101.9 ks 


1.52 ± 0.14 


(4.8 ± 0.4) x lO -11 


12-60 keV 


ISGRI 


24 March 2003-19 August 2009 


508 ks 


1.90 ± 0.20 


(1.2 ± 0.2) x 10-'° 


20-200 keV 


BAT 


December 2004-May 2009 


54 months" 


1.99 ±0.09 


(1.05 ± 0.05) x 10" 10 


15-150 keV 


LAT 


4 August 2008-20 April 201 1 


83.7 Ms 


2.4 ± 0.2 


(6 ± 2) x 10~ 12 


> 100 MeV 



Notes. (a) Elapsed time 



contamination by starburst emission. 

Using our best fit of the X-ray spectrum, we were able to con- 
strain some of the parameters. The indices of the power-law 
are tied to the power-law index of the X-ray spectrum: the first 
power-law index is p\ = 2 (below E^). Starting with this pa- 
rameter and the initial parameters of the code, we optimised the 
model. We focused on the energies that characterise the electron 
population to define the shape of the peaks, and to change the 
overall energy output in the SED we adjusted the radius of the 
emitting region and its magnetic field, and the Doppler factor. 
We found that the break in the electron power-law is insignificant 
and we therefore used a single power-law with index p\ - p2- 
Furthermore, we found the minimum energy of the electron dis- 
tribution to be E m i n = 5.6 x 10 6 eV and the maximum energy 
£max = 6.8 x 10 9 eV. The Doppler factor is S — 14, the radius 
of the spherical emission-volume is 7? s = 2 x 10 16 cm, and the 
magnetic field is B — 4 x 10~ 4 G. The parameters that we used 
to fit the data can be found in Table [2] and the resulting plot 
is presented in Figure [3] The addition of an external Compton 
component does not improve the model fit. 



/ 



1e+18 1e+20 1e+22 1e+24 



Fig. 3. SED showing 3C 111 unabsorbed fluxes and the SSC 
one-zone model. Crosses indicate data extracted and analysed 
in this work. We have added archival and deabsorbed IR and ra- 
dio points from NED. The line shows the SSC model. See Table 
[2]for the values used. 



4. Discussion 

To evaluate the physical processes that dominate the high-energy 
emission from X-ray to y-rays in 3C 1 1 1, we studied the X-ray 
spectrum and the broad-band SED. 

4.1. X-ray spectrum 

We constrained the high-energy cut-off of 3C 111 by ap- 
plying the pexrav model, which describes an exponentially 
cut-off power-law spec t rum r eflected from neutral material 
(Magdziarz & Zdzi arskil 1 1995b in xspec We found a cut- 
off E cut = 227^f k e v using data from INTEGRAL/1SGRL, 
SuzakufXlS and PIN, and Swift/BAI. An indication of the cut- 
off can already be seen in the power-law indices of the individ- 
ual fits, which go from r ~ 1 .6 in the soft X-rays to ~ 1 .9 in 
the hard band (Table [TJ. We also applied the physical compPS 
model to the same data, which yielded an electron temperature 
kT e = 91^keV that can be related to the cut-off energy via 
E cut 3kT e . The electron temperature that we found using the 
compPS model is consistent with the cut-off energy that we de- 
rived using the pexrav model. 

The value we found for the high-energ y cut-o ff is co nsis- 
tent with t he u pper limits reported by Dadina (2007 ) an d 
Ballo etalJteOllh and the measurement by Moli na et alT (l2009). 
Rivers et all (1201 ll) found, using data from RXTE from 3 keV to 
> 100 keV, that the addition of a cut-off does not improve their 
fit of 3C 111. This might be because we used a combination of 
different instruments and a wider energy range for the spectrum 
presented here. 

A high-energy cut-off is a typical property of the high-energy 
spectra of Seyfert galaxies, and the value we derived for 3C 111 
is in the expected range for cut-off en ergies observed in Seyferts 
(Beckmann et al . 2009; Dadina 2008). 

We measured the reflection in 3C 111 using both the compPS 
and the pexrav model to be R = 1.8^ and R = 0.7 + 0.3, 
respectively. Since these models are very different (pexrav is 
a descriptive model and compPS is a physical model), we ex- 
pect to find different measures of the reflection. iBeckm ann et al. 
(1201 ll) applied both models to the hard X-ray spectrum of the 
radio galaxy Cen A and found that the reflection in the case of 
compPS is also slightly higher though consistent with the reflec- 
tion found using the pexrav model. 

The reflection scaling factor we derived, R = 0.7 + 0.3, using the 
pexrav model can be compared to earlier work using the same 
m odel. Ou r result s are consistent with an upper limit reported 
bv iDadinal d2007l). as w ell as with the con straints obtained by 
Moli na et al] d2009l) and iBallo et al.1 (l201ll) . The latter also re- 
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Table 2. Parameters used to model the SEP of 3C 111 using th e code bv lKrawczvnski e t al. (2004). F or comparison, we in clude the 
fit parameters for the core of Centaurus A (lAbdo et al.ll2010bh and Markarian 421 (low flux state, Blazeiow ski et al.l l2005). 



Symbol 


Description 


3C 111 


Cen A 


Mrk421 


6 


Doppler factor 


14 


1.0 


10.0 


R s 


Radius of the emission volume [cm] 


2 x 10 16 


3.0 x 10 15 


7.0 x 10 15 


B 


Magnetic field of emission volume [G] 


4 x 10~ 2 


6.2 


0.405 


F 


Minimum energy of the electron distribution [eV] 


5.6 x 10 6 


1.5 x 10 s 


3.2 x 10 6 


Envdn. 


Maximum energy of the electron distribution [eV] 


6.8 x 10 9 


5.0 x 10 13 


1.7 x 10" 


Ebr 


Break energy of the electron distribution [eV] 




4.0 x 10 s 


2.2 x 10 10 


Pi 


Spectral index of electron spectrum (Emm to E^) 


2 


1.8 


2.05 


Pi 


Spectral index of electron spectrum (E bl to E m . ix ) 




4.3 


3.6 


logL 


Bolometric luminosity (3-200 keV) [ergs -1 ] 


44.7 


42.5" 


45.5* 



Notes. (fl) Beckmann et al. (2011); {h) Lichti et al. (2008) 



ported a significantly lower result, R = 0.19^^, when using 
only soft X-ray data between 0.4 keV and 10 keV. This is to be 
expected since the reflection compo nent depends on the energy 
band used to observe it. In addition R ivers et al. I d2oTiT) did not 
detect a reflection component when applying the pexrav model 
between 3 keV and > 100 keV. They constrained the model so 
as to avoid a high-energy cut-off, which might explain the differ- 
ence from our result. 

The value we measured is also consistent with the average 
Seyfert properties values of R = 1.2*^ for Se yfert 1 and 
R = l.l^ol for Seyfert 2 galaxies (iBeckmann et al.ll2009l) . 
For the iron line at 6.11 keV, we foun d that its EW = 85 ± 
1 1 eV, which is similar to the results of IB alio et al.l (120111) and 
Tom besi et al.l (l2010h . Smaller values between EW > 33 eV and 
EW > 40 eV were reported bv lTombesi et al.l (1201 ll) . Since the 
source is variable, it is likely that the iron line EW varies de- 
pending on the continuum. Figure [4] shows the correlation be- 
tween the flux between 4 keV and 10 keV and the EW of the 
iron line for several observations. The correlation between the 
flux and the EW has a significance of > 95%, a Pearson's test 
coefficient of r = -0.3, and a Spearman's rank correlation co- 
efficient r s = -0.7, showing a believable correlation between 
the flux level and the EW. When we did not consi der the large 
EW m easurement based on RXTE data presented bv lRivers et al.l 
(1201 ll) . w e get a highly significant cor relation with a probability 
of > 99%. To mbesi et al.l(l2010il201 lb also report the power-law 
indices changing for the different data sets. Similarly, using a 
different data set. lBallo et all d201 ll) derived an EW = 38^' eV. 
Since there is a difference of several months between both data 
sets this value can again be explained by the source variability. 

We found th at the EW is smaller than the average value for 
Seyfert galaxies iDadinal (120081) . at EW = 448 ± 67 eV, which 
can be compared to the average for Seyfert 1 galaxies of EW = 
222 + 33 eV and for Seyfert 2 galaxies of EW = 693 + 195 eV. 
The correlation between EW of the iron line and the underlying 
power-law flux indicates that the line flux does not vary signif- 
icantly, whereas the continuum varies by a factor of ~ 5. This 
can be naturally explained if the majority of the continuum flux 
is produced in a non-thermal and variable process. The small 
value of the EW of the iron line may indicate that the X-ray 
spectrum of 3C 111 does not only contain a thermal, but also a 
non-thermal contribution. 

We also fitted the spectrum with the physical compPS model, de- 
riving a temperature for the electron cloud kT e — 91+?? ke V, a 



I 



Flux (101 1 ergs/cm /s) 

Fig. 4. Flux of the source versus the EW of the iron line . 
Data taken fromlTombesi et al.l (I2010U201 lh . lBallo et all d201 ll) 
Rive rs et all (|201 ll) . Lewis et all d2005l) . and [Eracleou s et al.l 
(2000) and this work are plotted as points with error bars. We 
used a linear regression fit to derive a correlation between the 
two parameters. The fit is plotted with grey dashed lines indicat- 
ing the lcr error region. 



Compton parameter y = 0.6 + 0.1 (which corresponds to an op- 
tical depth of r = 0.8), and a reflection component/? — 1-8^. 
NGC 4151 is a well-s t udied case of a Seyfert core at hard X- 
rays. IBeckmann et al. I J2005I) fitted the compPS model to the 
spectrum of the Seyfert galaxy NGC 4151 across the range 
2-300 keV. They found for the electron plasma temperature 
kT e = 94+{ keV, a reflection component R = 0.72 + 0.14 
and an optical depth of r = 1.3+gH. The hard X-ray spec- 
trum (3-1000 keV) of Centaurus A, which is another nearby 
y-ray detected radio ga laxy, was also modelled with compPS by 
IBeckmann et al.l (1201 ll) . This yielded an electron plasma temper- 



ature of kT e = 206 + 62 keV, a Compton parameter y = 0.42^ ^ 
(which corresponds to an optical depth of r = 0.26), and a reflec- 
tion component of R = 0.12+° ^. NGC 4151 is optically thick, 
whereas both Cen A and 3C 111 are not. The optical thick- 
ness of NGC 4151 depends on the flux state and can d ecrease 
to t =~ 0.3 - 0.6 in the dim state (Lubins ki et al.ll20Tol) . when 
the electron temperature can increase to T e ~ 180 - 230 keV. 
The optical depth of 3C 111 is similar to that of the dim state 
of NGC 4151, but the electron temperature is closer to the elec- 
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tron temperature of the brighter state. The reflection measured in 
3C 111 is higher than th ose found in both NGC 4151 and Cen A. 
IChatteriee et al.l d201 lb concluded that the X-ray spectrum of 
3C 111 is of thermal inverse Compton origin, based on a cor- 
relation between the optical and X-ray fluxes, as well as a weak 
optical polarisation and a smaller variance in the optical than X- 
ray flux on shorter timescales. They concluded that these find- 
ings are consistent with a reprocessing model where the X-rays 
are mostly produced by inverse Compton scattering of thermal 
optical/UV seed photons from the accretion disc. 
The cut-off we measured using the pexrav model may have 
either a thermal or a non-thermal origin. It can be interpreted 
as a high-energy cut-off measured in Seyfert galaxies. In the 
case of 3C 111, the cut-off value is also comparable to the typ- 
ical cut-offs found in Seyfert galaxies. The cut-off observed in 
the pexrav model can also be caused by non-thermal inverse- 
Compton scattering processes. In this case, the smooth turn-over 
of the inverse Compton branch towards higher energies appears 
as a cut-off in the pexrav model. While in the thermal inverse- 
Compton case the cut-off is exponential, the non-thermal spec- 
trum is curved in such a way that y-ray and VHE emission can- 
not be excluded. It is impossible to differentiate between an ex- 
ponential and a simple cut-off in the hard X-rays based on the 
data at hand, thus both the non-thermal and the thermal interpre- 
tations of the curved X-ray spectrum are still valid. The reflec- 
tion component is a property of Seyfert galaxies and a result of 
thermal processes. The value of the one we detect is in the typ- 
ical range for reflection in Seyfert galaxies. The EW of the ob- 
served iron line is smaller than expected for Seyfert galaxies and 
is also variable, implying that the continuum is similarly vari- 
able. We therefore suggest that, while the thermal processes in 
the X-ray band in 3C 111 seem to dominate, we cannot exclude 
a non-thermal contribution. 



4.2. y-rays 

3C 111 was proposed as a counterpart of the y-ray source 
3EG J0416+3650 in the third CGRO/EGKET c atalogue, even 
thoug h it fell outside the 99% pro bability region dHartman et alj 
1999). A re-analysis of the data (Sgu era et al1l2005h concluded 
that 3EG J0416+3650 is likely associated with 3C 111. There 
are no other plausible counterparts in the EGRET error re- 
gion, which is also larger than previously thought because the 
quoted errors were statistical only and did not take into ac- 
count the larger systematic errors caused by inaccuracies in 
the Galactic diff use model. The EGRET data were re-analysed 
dHartman et al.ll2008l) when it was found that 3 EG J0416+3650 
is likely the result of the blending of more than one source. One 
of these components (detected only above 1 GeV) can be asso- 
ciated with 3C 111. Furthermore , 3C 111 was incl uded in the 
first Fermi/LAT source catalogue ( Ab do et a ni20T0ah with a sig- 
nificance of 4.3cr. In the second Fermi/LAT source catalogue, 
3C 111 was excl uded because the source was no longer signif- 
icantly detected dNolan et al.l 12012 ). Howeve r, 3C 111 is very 
likely to be variable ( Ackermann et alJl201 lh and therefore no 
longer detectable in the second year. This was also confirmed by 
Kataoka et al.l d201 ll) . who found a signific ance >5<x for 3 C 111 
using 24 months of Fermi/LAT data and by Grandi et all d2012l) 
who found that the y-ray emission is not persistent, but flaring 
and associated with the ejection of bright radio knots. 
Since 3C 1 1 1 was detected by CGRO/EGRET and Fermi/LAT 
at different epochs, we assumed that the source is a variable y- 
ray emitter. We included data from the first Fermi/LAT catalogue 
where the source was detected significantly. Our analysis gives a 



power-law index similar to that given in th e first catalogue with 
a comparable flux level (Abdo et al. 20 10ah . 



4.3. Spectral energy distribution 

To construct the SED of a variable source, simultaneously ac- 
quired data are necessary in order to ensure that we make an 
unbiased measurement. The data used in the SED of 3C 111 
were not acquired simultaneously, which can affect the results 
in the sense that the SED is composed of measurements from 
different spectral states. This effect is more severe for very vari- 
able objects, whereas at most waveb ands 3C 111 is moder- 
ately variable. We refer, for example to lBeckmann et all |2007) 
who found no significant variability in hard X-rays for their 
study of 9 months of S wift/BAT data. For the radio domain, 
iGrossberger et alJ d2012) present decade-long light curves show- 
ing variations of up to a factor of 2 which we can consider as in- 
significant variability in the context of the SED. For other radio 
galaxies, it has also been shown that the use of non-simultaneous 
data to construct the SED does not introduce any significant bias . 
For example, the radio galaxies NGC 1275 dAnton et alj|2004h 
and Pictor A dBrow n & Adams 2012) were also similarly anal- 
ysed using time-averaged SEDs. The results for the physical pa- 
rameters should thus not be significantly affected by the moder- 
ate variability of 3C 111. 

The archival radio data that we used to model the SED were not 
well-represented by our SSC model. For observatories that oper- 
ate at different wavebands, there are large differences in the data 
acquired in terms of the field of view and resolution. The resolu- 
tion depends on the wavelength via RL = A/D, where RL is the 
resolution, A the wavelength, and D the dish diameter. It is there- 
fore possible that we probe different regimes in the broad-band 
SED, something that the one-zone model does not account for. 
In SSC models, the relativistic electron-energy distribution is of- 
ten assumed to be a broken power-law with a break energy £br, 
a power-law index before the break of p\ < 3 and a power-law 
index after the break of P2 > 3 (see for example [ Ghisellin feTaT] 
ll996r.lTavecchio et alJl998l:lKrawczvnski et alJ2004h . The elec- 
tron population is confined to region of radius R s and a magnetic 
field B that moves with a Doppler factor 6 along the jet. The elec- 
trons emit synchrotron radiation and inverse Compton radiation 
using the synchrotron photons as seed photons. For electrons 
with energies below the break energy E\, T , cooling of the elec- 
trons is the dominant emission source and in the regime after the 
break, the escape of electrons from the source dominates. The 
peak synchrotron power is emitted by electrons at the break en- 
ergy dTavecchio et al.ll998h . In the case of 3C 1 1 1, we found that 
using a broken power-law gave a break energy of E\, r =1.6 GeV 
and power-law indices of p\ -2 and p2 = 2.2. The break be- 
tween the two power-laws is insignificant. We therefore applied 
a single power-law with an index p = 2, based on the power-law 
index of the X-ray spectrum. In the case of a single power-law, 
the maximum synchrotron power is emitted near the maximum 
energy of £ max =* 7 GeV. 

The energies that define the electron distribution, E m [ n and E m!lx 
influence the shape of the synchrotron and inverse Compton 
peak and were chosen by ourselves empirically. If the mini- 
mum energy, E m [ n , decreases, the maximum of the peak occurs 
at lower frequency and the slopes are less steep. Increasing the 
E m [ n increases the depth between the synchrotron and inverse 
Compton peak and steepens the rise and fall of the peaks. The 
maximum energy, E max , is tied to the peak frequencies of both 
the synchrotron and inverse Compton peak: increasing this value 
would increase the peak frequencies. 
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lAbdo et al.l (l2010bl) modelled the SED of the core of the ra- 
dio galaxy Cen A, using simultaneous data from the radio up 
to the y-ray regime, applying a single-zone SSC. Comparing 
the parameters for this model and our own (see Table |2j, it is 
clear that the models have significantly different parameter val- 
ues. The magnetic field has a much higher value in the case 
of Cen A (B = 6.2 G), than the value we found for 3C 111 
(B = 0.04 G). The higher magnetic field increases the result- 
ing flux because the synchrotron power depends on the mag- 
netic field. In contrast, the Doppler factor used to model Cen A 
is low (5 = 1) compared to the value we found for 3C 111, 
6 = 14. A lower Doppler factor means that the emission is less 
boosted and therefore appears less energetic. There are also three 
orders of magnitude difference between the emission volumes: 
for Cen A, R s — 3 x 10 15 cm was used, whereas we assumed 
R s = 2 x 10 16 cm for 3C 111. Since this defines the amount of 
particles this also causes the flux to be lower for Cen A than for 
3C 111. The higher value used for the magnetic field is attenu- 
ated by the lower Doppler factor and radius for Cen A, which 
causes the overall output to be lower, as expected for the less lu- 
minous source. 

Pictor A is another FR-II radio galaxy tha t was detected by 
Fermi/LAI as reported by iBrown & Adamsl d2012l) . who mod- 
elled the SED of one of the hot-spots in the radio-lobe of Pic A 
with a one-zone SSC model using the same code that we applied. 
They concluded that it is impossible to describe both the X-ray 
and y-ray emission using the SSC model. Since the y-ray emis- 
sion is very variable it is likely to originate from the jet, whereas 
the X-rays originate from the hot-spot. This appears to indicate 
that more than one zone i s needed to model the entire SED. 
Blaz eiowski et al.l (120051) also used the Rrawczynski code to 
model the high-energy peaked BL Lac object Mrk 421, assum- 
ing only an SSC component. The parameter values they derived 
are shown in Table [2] Both the Doppler factor (6 = 10) and the 
emitting region radius (R s — 7 x 10 15 cm) chosen for Mrk 421 
are smaller than those chosen for 3C 111, but are larger than 
those of Cen A. Similar to Cen A, the magnetic field value for 
Mrk 421, B - 0.4 G, is a factor of 10 stronger than for 3C 111. 
The smaller radius and Doppler factor decrease the flux, but ow- 
ing to the strong magnetic field the flux in the SED for Mrk 421 
is higher than the flux of 3C 111. 

We conclude that the overall emission from 3C 111 can be mod- 
elled with a simple synchrotron self-Compton model, where no 
additional thermal Compton component is needed. Since the X- 
rays appear to be (mostly) of thermal origin, the SED provides 
an upper limit to the non-thermal emission in the X-ray band. 
There may be a thermal component in the SED, but with the 
current data we are unable to disentangle the possible thermal 
component from the overall non-thermal emission. 

5. Conclusion 

The origin of th e high-energy emission from non-blazar AGNs 
remains unclear. iMarscher et al.l {2002) suggested that with the 
acceleration of the inner regions of the accretion disc a shock 
front will stream along the je t and the expulsion of a superlu- 
minal bright knot will follow. iGrandi et al.l (1201 2|) show, using 
Fermi/LAiT data, that the GeV emission of 3C 111 appears to 
originate from a compact knot confined to within 0. 1 pc. This 
knot is clearly separate from the core and placed 0.3 pc from the 
central engine. 

By analysing the X-ray spectrum and the broad-band SED, we 
have studied the nature of the high-energy emission of the ra- 
dio galaxy 3C 111. We have presented an X-ray spectrum be- 



tween 0.4 keV and 200 keV using data of 3C 111 acquired by 
several instruments and showed that the best-fit model is an ab- 
sorbed cut-off power-law with both a reflection component and a 
Gaussian component to account for the iron line at 6.4 keV. The 
values we found for the reflection and high-energy cut-off are 
similar to those found in Seyfert galaxies, which would indicate 
that there is a thermal core visible. The cut-off can also orig- 
inate from non-thermal processes and the EW of the iron line 
is variable and smaller than expected for Seyfert galaxies. We 
therefore conclude that the X-ray spectrum is mainly of thermal 
origin, but there may be a small non-thermal contribution. 
Using the X-ray spectrum, together with y-ray data from 
FermifLPH and archival deabsorbed radio and infrared data, we 
modelled the broad-band SED of 3C 111 using a single-zone 
synchrotron self-Compton model. This model is non-thermal 
and we also did not need to include an additional thermal com- 
ponent to model the SED. Since the X-ray emission is likely to 
have a combined thermal and non-thermal origin, the SSC model 
we used may overestimate the non-thermal contribution in the X- 
ray band and should therefore be considered an upper limit. 
In conclusion, it seems that the high-energy emission from 
3C 111 consists of both thermal and non-thermal components. 
In the X-ray spectrum, the thermal components manifest them- 
selves in terms of an iron line and reflection. The non-thermal 
component is visible through the variability in the EW of the 
iron line. The high-energy cut-off can be the result of either ther- 
mal or non-thermal inverse Compton scattering, but our present 
spectrum does not allow us to distinguish which process is occur- 
ring. The broadband SED can be modelled with a non-thermal 
model, but it is possible there is a thermal component that we are 
unable to discern with the current data set. 
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